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Abstract 

The  end-to-end  performance  achieved  by  an  adaptive  optical  (AO)  imaging  sys¬ 
tem  is  determined  by  a  combination  of  the  residual  time-varying  phase  distortions 
associated  with  atmospheric  turbulence,  and  the  quasi-static  unsensed  and  uncor- 
rectable  aberrations  in  the  optical  system  itself.  Although  the  effects  of  these  two 
errors  on  the  time  averaged  Strehl  ratio  and  the  time  averaged  optical  transfer  func¬ 
tion  (OTF)  of  the  AO  system  are  not  formally  separable,  such  an  approximation  is 
found  to  be  accurate  to  within  a  few  per  cent  for  a  range  of  representative  resid¬ 
ual  wave-front  errors.  In  these  calculations,  we  have  combined  static  optical  system 
aberrations  and  time-varying  residual  phase  distortion  characteristic  of  deformable 
mirror  (DM)  fitting  error,  wave-front  sensor  (WFS)  noise,  and  anisoplanatism.  The 
static  aberrations  consist  of  focus  errors  of  varying  magnitudes  as  well  as  a  combina¬ 
tion  of  unsensed  and  uncorrectable  mirror  figure  errors  derived  from  modeling  by  the 
Gemini  8-meter  Telescopes  project.  The  overall  Strehl  ratios  and  OTF's  due  to  the 
combined  effect  of  these  error  sources  are  well  approximated  as  products  of  separate 
factors  for  the  static  and  time-varying  aberrations,  as  long  as  the  overall  Strehl  ratio 
due  to  both  errors  is  greater  than  about  0.1.  For  lower  Strehl  ratios  the  products 
provide  lower  bounds  on  the  actual  values  of  the  Strehl  ratio  and  OTF.  The  speckle 
transfer  function  (STF)  is  also  well  approximated  by  a  product  of  two  functions,  but 
only  where  AO  compensation  is  sufficiently  good  that  speckle  imaging  techniques  are 
usually  not  required. 

OCIS  codes:  010.1080,  010.4850.  100.1830, 110.6150, 110.6770 


1.  Introduction 

The  existing  literature  suggests  that  most  adaptive  optics  (AO)  performance  analysis  for 
large,  ground  based  telescopes  is  partitioned  into  separate  treatments  of  the  phase  distor¬ 
tions,  or  aberrations,  introduced  by  atmospheric  turbulence,  and  aberrations  introduced  by 
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the  telescope  optics  themselves. Although  much  expense  and  modeling  effort  is  applied 
to  characterize  these  two  wave-front  error  sources  separately,  overall  optical  system  perfor¬ 
mance  is  frequently  estimatedJ)y  simply  multiplying  the  corresponding  Strehl  ratios.  For 
A0-c6mpensated  wave-fronts  with  residual  errors,  this  multiplicative  approximation  is  not 
formally  exact:  since  the  costs  and  performance  specifications  associated  with  large  tele¬ 
scopes  are  so  challenging,  it  may  be  prudent  to  model  the  interaction  of  turbulence- induced 
phase  distortions  and  telescope  aberrations  with  greater  precision.  This  paper  summarizes 
the  results  of  such  a  study.  For  a  range  of  representative  atmospheric  phase  distortions 
partially  compensated  by  adaptive  optics  and  several  sample  telescope  aberrations,  the  mul¬ 
tiplicative  estimate  of  the  overall  system  Strehl  ratio  is  a  remarkably  good  approximation 
as  long  as  this  quantity  is  larger  than  about  0.1.  At  somewhat  lower  performance  levels 
the  multiplicative  approximation  can  be  used  as  a  moderately  conservative  lower  bound  on 
the  actual  Strehl  ratio.  Similar  conclusions  apply  to  approximating  an  “end-to-end”  optical 
transfer  function  (OTF)  which  accounts  for  all  of  the  wave-front  errors  in  the  system  as 
a  product  of  two  terms  computed  separately  for  the  effects  of  the  atmospheric  turbulence 
and  the  telescope  optics.  For  higher-order  moments  of  the  OTF  such  as  the  speckle  trans¬ 
fer  function  (STF)  the  multiplicative  approximation  begins  to  fail  at  a  somewhat  higher 
Strehl  ratio,  but  reasonable  accuracy  can  still  be  achieved  at  Strehl  ratios  of  about  0.3  if  the 
residual  wave-front  aberrations  are  distributed  over  a  range  of  spatial  frequencies. 

The  validity  of  these  approximations  simplifies  the  design  and  use  of  AO  systems  for 
large  telescopes  in  several  ways.  During  system  design  an  overall  system  OTF  specification 
may  be  factored  into  two  allocations  for  the  telescope  optics  and  the  AO  system  which 
may  be  addressed  separately.  When  AO-compensated  imagery  is  to  be  postprocessed  using 
deconvolution  of  the  OTF,  the  end-to-end  system  OTF  can  be  approximated  as  a  product 
of  (a)  the  atmospheric  OTF  as  estimated  from  closed-loop  wave-front  sensor  (WFS)  data^ 
and  (b)  the  optical  system  OTF  as  calibrated  using  a  local  reference  source.  If  speckle 
imaging  techniques'*  are  used,  the  system  STF  can  also  be  approximated  by  the  product  of 
an  atmospheric  STF  (estimated  from  WFS  data)  and  an  optical  system  STF. 

The  phase-only  Fourier  optics  model  used  for  this  study  is  described  in  section  2.  The 
end-to-end  optical  system  OTF  is  proportional  to  the  autocorrelation  of  the  complex-valued 
function  U.  the  optical  field  obtained  in  the  exit  pupil  of  the  telescope  from  a  monochromatic 
point  source  after  propagation  through  the  atmosphere,  the  AO  system,  and  the  telescope 
optics.  The  pupil  function  U  can  be  divided  into  three  parts:  The  real- valued  transmittance 
function  of  the  optical  system,  the  static  telescope  aberrations  <j>3  which  are  unsensed  or 
uncorrected  by  the  adaptive  optics,  and  the  time-varying  residual  atmospheric  distortions 
o,.  left  uncompensated  by  the  AO  system.  The  time-averaged  OTF  of  the  system  is  obtained 
by  averaging  over  the  statistics  of  the  atmospheric  phase  profile  (f)^.  which  is  assumed  to 
be  a  zero-mean,  normally  distributed  random  variable.  The  time  averaged  Strehl  ratio  S 
can  be  computed  as  the  integral  of  the  OTF.  Without  adaptive  optics,  the  stationarity  of 
turbulence-induced  phase  distortions  allows  the  OTF  to  be  written  as  a  product  of  two 
terms  which  depend  separately  upon  and  the  statistics  of  (j)^.  For  atmospheric  turbulence 
which  has  been  partially  compensated  by  AO  stationarity  may  not  be  assumed,  however, 
and  such  a  decoupling  is  not  mathematically  valid.  Experience  and  intuition  still  suggest 
that  for  small  phase  errors  the  combined  effect  of  and  should  be  well  approximated  by 
the  product  of  separately  computed  transfer  functions  or  Strehl  ratios.  For  speckle  imaging. 
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however,  the  form  of  the  STF  casts  doubt  on  the  validity  of  a  product  approximation  for 
the  end-to-end  transfer  function  even  without  AO  compensation. 

Sample  results  quantifying  the  accuracy  of  these  approximations  are  summarized  in  sec¬ 
tion  3  for  a  range  of  representative  telescope  optical  aberrations  and  residual  turbulence- 
induced  phase  errors.  Three  varieties  of  residual  turbulence-induced  errors  driven  primarily 
by  deformable  mirror  (DM)  fitting  error,  wave-front  sensor  (WFS)  noise,  and  anisoplanatism 
are  combined  with  three  different  magnitudes  of  a  static  focus  error  for  a  sample  Shack- 
Hartmann-based  AO  system  with  8x8  WFS  subapertures  and  9x9  DM  actuators.  Two 
more  specific  cases  are  also  considered  using  atmospheric,  AO,  and  aberration  parameters 
corresponding  to  the  Gemini-North  8-meter  telescope.®’®  The  results  of  these  calculations 
form  the  basis  for  the  conclusions  given  in  the  first  paragraph  above. 

Section  4  presents  sample  sky  coverage  results  obtained  using  the  Gemini  telescope  OTF’s 
as  computed  in  section  3.  Here  “sky  coverage”  is  defined  as  the  function  P{S  >  Sq).  the 
fraction  of  a  specified  region  of  the  sky  for  which  the  telescope  and  AO  system  achieve 
an  end-to-end  performance  level  5  at  least  equal  to  a  desired  threshold  Sample 

performance  metrics  for  sky  coverage  include  time  average  Strehl  ratio  and  the  coupling 
efficiency  for  feeding  an  optical  fiber  or  spectrometer  slit.  This  function  is  determined  by 
a  combination  of  guide  star,  atmospheric,  science  instrument,  AO  system,  and  telescope 
parameters,  including  the  unsensed  and  uncorrectable  aberrations  in  the  telescope  optics. 
For  sample  Gemini-North  parameters  we  find  that  the  reduction  in  sky  coverage  due  to 
telescope  aberrations  for  a  specific  instrument,  such  as  a  spectrometer,  can  depend  upon 
the  shape  of  the  aberrated  point  spread  function  (PSF)  as  well  as  upon  its  Strehl  ratio. 
In  this  case  the  end-to-end  OTF  of  the  AO  system,  accounting  for  both  atmospheric  and 
telescope  error  sources,  must  be  used  to  compute  accurate  sky  coverage  estimates.  Sky 
coverage  estimates  derived  by  simpler  scaling  laws  can  be  significantly  pessimistic. 

Finally,  section  5  is  a  brief  summary  of  the  principal  results  of  the  paper. 


2.  Analysis  Summary 


Using  Fourier  optics  methods,^®’^'  the  monochomatic,  incoherent,  isoplanatic  imaging  per¬ 
formance  of  an  optical  system  can  be  modeled  as  a  spatial  filtering  operation  by  the  transfer 
function  OTF(k).  This  tranfer  function  is  described  by  the  equation 


OTF(/c) 


fdr  Uir)U‘{r-\K) 
fdr\U{r)\^ 


(2.1) 


Where  k  is  a  two-dimensional  spatial  frequency  variable,  A  is  the  wavelength  of  the  light, 
the  integration  variable  r  denotes  coordinates  in  the  telescope  exit  pupil,  and  the  function 
U  is  the  optical  field  obtained  in  the  exit  pupil  when  imaging  a  point  source.  If  atmospheric 
scintillation  is  neglected,  the  field  U  is  related  to  the  residual  turbulence-induced  phase 
distortion  <p^.  left  uncorrected  by  the  AO  system  by  the  formula 


U(r)  =  P(r)exp[f?i„(r)] ,  (2.2) 

where  P  is  the  field  in  the  exit  pupil  with  no  atmospheric  turbulence.  The  notation  “<l>v^ 
has  been  selected  to  suggest  that  the  residual  turbulence-induced  phase  distortion  is  a  time- 
varying  quantity,  and  this  function  is  assumed  to  be  a  zero-mean,  normally  distributed 
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random  variable.  The  turbulence- free  field  P  is  referred  to  as  the  generalized  pupil  function 
of  the  telescope  and  may  be  written  in  the  form 

P{r)  =  Poir)exp[iUr)]:  (2.3) 


where  Pq  is  a  real-valued  optical  transmittance  function  and  <f>a  is  the  static  or  quasi-static 
wave-front  error  due  to  uncorrectable  or  unsensed  aberrations  in  the  telescope  optics.  Ex¬ 
amples  of  such  errors  include  errors  at  unobservably  high  spatial  frequencies  on  the  primary 
mirror,  and  aberrations  in  optical  elements  within  instruments  which  are  not  included  in 
the  WFS  optical  path.  Substituting  Eq.'s  (2.2)  and  (2.3)  back  into  Eq.  (2.1)  yields  the 
relationship 


OTF(k;)  = 


fdr  Po{r)Po(r  -  AK)exp  {f  [0*(r)  -  <j>s{r  -  AK)]}exp{i  [<;6„(r)  -  <j>y{r  -  A#c)]} 


f dr  P^{r) 


(2.4) 


between  the  optical  system  OTF  and  the  wave-front  aberrations  <l>s  and  induced  by  the 
telescope  optics  and  uncompensated  atmospheric  turbulence. 

For  long  exposure  imaging  it  is  the  time  averaged  value  (OTF(k))  which  is  of  interest, 
where  the  notation  (•  •  •)  denotes  ensemble  averaging  over  the  statistics  of  atmospheric  tur¬ 
bulence  and  WFS  measurement  noise.  Since  the  telescope  aberration  function  is  assumed 
to  be  static,  only  the  final  term  in  the  numerator  of  Eq.  (2.4)  is  a  random  quantity.  Be¬ 
cause  the  uncorrected  turbulence-induced  aberration  (j)^  is  modeled  as  a  zero-mean,  normally 
distributed  random  variable,  the  expected  value  of  this  term  takes  the  form 

(e.xp  {?■  [d>„(r)  -  4>„{r  -  A/c)]})  =  exp  ([<;i„(r)  -  -  Ak)]^)} 

=  exp  [— 5l)(r,  r  —  A/c)] .  (2.5) 

where  T>{r.r  —  Ak)  =  ^[^ir(**)  —  is  the  structure  function  of  the  residual 

turbulence-induced  phase  error  </»„.  The  time  averaged  OTF  of  the  optical  system  is  therefore 
given  by  the  expression 


(OTF(k)) 


fdr  Po{r)Po{r  —  AK)exp  {» [<?!>5(r)  —  <f>s{r  —  Ak)]}  exp  —\V{r,r  —  AK;)j 

idrPiir) 


(2.6) 


The  point  spread  function  (PSF)  of  the  optical  system  and  related  quantities  can  be  now 
evaluated  in  terms  of  the  time  averaged  OTF.  In  particular,  the  time  averaged  Strehl  ratio 
5  of  the  optical  system  can  be  computed  using  the  formula 


5  =  J  dK  (OTF(k))  .  (2.7) 

From  the  last  two  equations  it  is  clear  that  the  time  averaged  OTF  and  the  Strehl 
ratio  5  are  functions  of  the  static  telescope  aberration  and  T>,  the  structure  function 
of  the  time-varying  phase  errors  (j)^.  It  is  convenient  to  abbreviate  this  relationship  using 
the  notation  (OTF(k))(^^  uj  and  For  general  values  of  (j)a  and  V  the  dependance  of 
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the  OTF  and  the  Strehl  ratio  upon  these  quantities  is  nonseparable.  except  in  the  special 
case  where  either  the  static  phase  difference  <l>sir)  —  <l>s{r  —  A/c)  or  the  structure  function 
T>{r.r  —  \k)  is  actually  a  function  of  Xk  alone.  This  is  indeed  the  case  for  stationary  models 
of  uncompensated  atmospheric  turbulence,  but  these  conditions  do  not  generally  hold  for 
the  residual  wave-front  errors  left  uncorrected  by  an  AO  system.  Experience  and  intuition 
suggest,  however,  that  multiplicative  approximation 

(2.8) 

should  still  be  reasonably  accurate  as  long  as  the  quantities  d**  and  T>  are  hot  too  large. 

For  the  Strehl  ratio  S.  the  corresponding  multiplicative  approximation  can  be  justified 
for  two  types  of  phase  errors.  First,  AO  systems  operating  under  benign  conditions  will 
leave  small  residual  phase  errors.  For  this  case  we  may  use^^ 

S  Js  1  -  (2.9) 


{OTF(k))|^„p,  a  {OTF(k)),„,„I 


(OTF(«))„.,„, 


where  cr^  is  the  variance  of  (j>  over  the  telescope  aperture.  While  this  approximation  would 
not  be  accurate  for  uncompensated  turbulence, the  assumption  of  small  residual  phase 
errors  justifies  its  use  in  this  case.  Substituting  (f>  =  <i>v  +  <f>s  into  Eq.  (2.9)  and  assuming 
that  the  two  terms  are  spatially  uncorrelated  yields 

ss  +  aj) 

ss  exp  [-(<7^  +  aj}] 

~  5(O.T>)‘5(0j.O);  (2.10) 

where  and  cr^  denote  the  variances  of  the  two  terms  of  the  wave-front  error.  For  sufficiently 
small  total  phase,  then,  the  overall  Strehl  ratio  may  be  approximated  by  the  product  of 
two  independent  terms,  regardless  of  the  spatial  distribution  of  the  phase  errors.  Secondly, 
consider  the  limit  case  of  residual  turbulence-induced  phase  errors  which  are  spatially  white, 
so  that  V{r.r  —  Xk)  =  2crl  for  all  r  and  all  nonzero  k.  Combining  this  assumption  with 
Eq.’s  (2.6)  and  (2.7)  yields 


5  =  J  dK  (OTF(k)) 


2  J<^’’-Po(»’)Po(r  -  AK))exp{i[?i>,(r)  -  <f>sir  -  Ak)]} 


(2.11) 


so  that  the  overall  Strehl  is  again  given  as  the  product  of  two  terms.  As  shown  by  the  results 
presented  in  Section  3  below,  one  of  these  two  assumptions  on  the  distribution  of  cr„  will 
generally  be  valid  for  an  AO  system  providing  any  useful  degree  of  atmospheric  turbulence 
compensation. 

Finally,  for  some  applications  higher  moments  of  the  OTF  may  be  of  interest  as  well. 
Low-order  AO  systems  designed  to  compensate  turbulence  with  very  dim  guide  star  refer¬ 
ences  can  fail,  in  severe  seeing  conditions,  to  provide  PSFs  in  the  so-called  “sharp  core” 
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regime,  where  there  is  a  reasonable  fraction  of  the  total  PSF  intensity  in  a  stable  diffraction- 
limited  central  core.  If  conditions  are  poor  enough,  speckle  imaging  techniques'^  can  be 
required  to  achieve  reasonable  stability  in  the  post-processed  image.^  While  many  deconvo¬ 
lution  methods  require  only  an  estimate  of  the  long-exposure  OTF,  speckle  imaging  without 
observation  of  a  calibration  star  would  require  WFS  estimates  of  higher  moments  of  the 
atmospheric  OTF,  such  as  the  speckle  transfer  function  (STF),  cross-spectrum,  and  bispec¬ 
trum.  For  this  reason,  we  also  consider  approximation  of  the  STF  by  the  product  of  the 
atmospheric  STF  and  the  optical  system  STF.  Under  the  same  conditions  used  to  develop 
the  OTF,  the  STF  is  given  by 

STF(k)  =  (|OTF(/c)|2) 

=  J  dr  J  dr'  Po{r)Po{r  -  XK)Po{r')Poir'  -  Xk) 

X  exp{-f[^,(r)  -  4>s{r  -  Xk)  -  <f>,{r')  +  <f>,{r'  -  A«)]} 

X  (exp{-t[«^„(r)  -  ^„(r  -  Ak)  -  +  ^^{r'  -  A/c)]}) .  (2.12) 


Because  the  expectation  is  over  phase  factors  which  are  functions  of  both  variables  of  integra¬ 
tion.  it  is  not  stationary  and  can  not  be  factored  outside  the  integrals.  If  AO  compensation 
is  good,  the  residual  atmospheric  phase  4>v  will  be  small,  but  in  this  case  the  resulting  images 
are  in  the  sharp  core  regime  where  we  assume  speckle  processing  would  not  be  used.  To  see 
whether  a  larger  domain  exists  over  which  the  STF  may  be  approximated  as  a  product  of 
factors,  we  will  consider  the  approximation 


STF(K),0,.t)) 


STF(k)(o,p) 


STF(k)(^„o) 

STF{k)(o.o) 


(2.13) 


The  extent  to  which  the  above  multiplicative  approximations  apply  for  a  range  of  repre¬ 
sentative  AO  imaging  scenarios  is  explored  in  the  following  section. 


3.  Sample  Numerical  Results 
A.  Cases  Considered 

The  formulas  developed  in  the  preceeding  section  have  been  applied  to  study  the  interac¬ 
tion  of  turbulence-induced  phase  errors  with  unsensed  and  uncorrectable  static  telescope 
aberrations  for  a  number  of  sample  AO  imaging  scenarios.  The  four  cases  considered  are 
summarized  in  Table  1,  The  first  three  cases  relate  to  a  generic  Shack-Hartmann-based  AO 
system  with  8x8  WFS  subapertures  and  9x9  DM  actuators,  with  residual  turbulence- 
induced  phase  errors  driven  by  DM  fitting  error,  WFS  noise,  and  anisoplanatism.  The 
techniques  used  to  evaluate,  the  structure  function  V  of  the  residual  turbulence-induced 
phase  distortions  for  each  of  these  cases  have  been  described  elsewhere.^®  The  static  tele¬ 
scope  aberration  is  a  simple  focus  error  for  all  three  of  these  cases;  in  practice  such  an  error 
could  result  from  either  an  optical  misalignment  or  by  miscalculating  the  effective  range  of 
a  laser  guide  star. 
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The  last  of  the  four  cases  is  based  upon  atmospheric.  AO.  and  telescope  parameters 
derived  from  design  studies  for  the  Gemini-North  8-meter  telescope.®'®  Here  the  turbulence- 
induced  residual  phase  errors  are  primarily  a  combination  of  fitting  error  and  focus  aniso- 
planatism  for  a  laser  guide  star  generated  in  the  mesospheric  sodium  layer.  An  estimate 
of  typical  unsensed  and  uncorrectable  optical  aberrations  for  Gemini-North  is  illustrated  in 
Fig.  1.  Uncorrectable  mirror  polishing  and  mounting  errors  have  been  studied  in  great  de¬ 
tail  for  Gemini,  while  the  noncommon  path  errors  between  the  WFS  and  science  instrument 
optics  remain  budgeted  allocations  at  this  time.  The  phase  profile  in  Fig.  1  is  a  sum  of 
measured  primary  mirror  mounting  and  polishing  errors,  a  very  small  simulated  mounting 
error  for  the  secondary  mirror,  and  a  “typical”  noncommon  path  error  represented  as  a  sum 
of  low-order  Zernike  modes.  This  phase  profile  yields  a  combined  telescope  and  instrument 
Strehl  ratios  of  about  0.70  at  a  wavelength  of  1.65 //m  .  and  a  more  benign  case  with  a 
Strehl  ratio  of  about  0.85  has  also  been  considered.  These  two  values  correspond  roughly  to 
the  Gemini  specification  for  the  noncommon  path  errors  and  the  smallest  errors  which  can 
realistically  be  expected. 


B.  Parametric  Results 

Fig.  2  illustrates  the  Strehl  ratios  computed  for  the  first  three  cases  described  in  Table  1. 
The  solid  curves  for  the  “exact  Strehl  ratio”  have  been  computed  using  Eq.  (2.7).  and  the 
dashed  curves  for  the  “multiplicative  approximation”  are  based  upon  Eq.  (2.10).  Results 
for  DM  fitting  error  combined  with  varying  amounts  of  static  focus  error  are  summarized  in 
Fig.  2a  for  subaperture  widths  equal  to  1,  2,  and  4  times  the  turbulence-induced  effective 
coherence  diameter  tq.  and  static  focus  errors  with  values  between  0  and  1  radian  RMS.  The 
multiplicative  approximation  to  the  overall  Strehl  ratio  has  a  worst-case  relative  error  of  no 
more  than  about  six  per  cent,  even  for  overall  Strehl  ratios  as  low  as  two  per  cent.  In  this 
last  case  the  approximation  is  accurate  even  with  relatively  large  residual  phase  distortions 
because  the  errors  are  high  spatial  frequency  and  approximately  white. 

.Analogous  Strehl  ratio  results  for  varying  amounts  of  random  phase  errors  due  to  WFS 
noise  and  anisoplanatism  are  plotted  in  Fig.'s  2b  and  2c.  Here  the  relative  errors  in  the 
multiplicative  approximation  remain  below  about  five  per  cent  for  Strehl  ratios  no  lower 
than  0.2.  and  no  more  than  about  ten  per  cent  for  Strehl  ratios  on  the  order  of  0.1.  These 
somewhat  larger  errors  in  the  approximation  occur  because  the  random  phase  errors  due  to 
WFS  noise  and  anisoplanatism  contain  nonnegligible  low  spatial  frequency  modes,  which  do 
not  necessarily  have  equal  power  everywhere  in  the  pupil. 

Results  on  the  accuracy  of  the  multiplicative  approximation  for  the  end-to-end  optical 
system  OTF  due  to  combinations  of  the  errors  considered  above  are  summarized  in  Fig.’s  3 
and  4.  The  curves  labeled  “exact  OTF”  and  “transfer  function  approximation”  have  been 
computed  using  Eq.’s  (2.6)  and  (2.8),  respectively.  These  plots  are  for  the  case  of  a  static 
focus  error  of  0.7  radians  RMS,  which  corresponds  to  a  Strehl  ratio  of  about  0.6  without 
atmospheric  turbulence.  As  illustrated  in  Fig.  3,  this  aberration  reduces  the  OTF  to  about 
45  per  cent  of  its  diffraction  limited  value  at  a  normalized  spatial  frequency  of  one-half  the 
diffraction-limited  cutoff  frequency.  The  OTF’s  computed  for  a  combination  of  DM  fitting 
error  and  this  static  focus  error  are  plotted  in  Fig.  4a.  The  agreement  between  the  actual 
OTF  and  the  multiplicative  approximation  is  excellent,  once  again  because  the  DM  fitting 
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error  has  so  little  low  spatial  frequency  content.  OTF  values  for  wave-front  errors  induced 
by  WFS  noise  or  anisoplanatism  are  summarized  in  Fig.'s  4b  and  4c.  Here  the  worst  case 
relative  error  between  the  actual  OTF  and  the  multiplicative  approximation  is  on  the  order 
of  *12  *to'  18  per  cent  for  cases* where  the  overall  end-to-end  Strehl  ratio  is  from  6  to  20  per 
cent,  and  about  five  percent  for  cases  where  the  end-to-end  Strehl  ratio  is  from  20  to  30 
per  cent.  This  level  of  agreement  should  be  adequate  for  developing  system  performance 
specifications  and  many  image  postprocessing  applications. 

Fig.  5  shows  STF’s  for  subaperture  widths  as  in  Fig.'s  2a  and  4a;  d  =  1,  2,  and  4 
times  To.  The  solid  curves  plot  the  exact  results  computed  using  Eq.  (2.12),  and  the  dashed 
curves  indicate  the  multiplicative  approximation  given  by  Eq.  (2.13).  Fig.  5a  plots  exact  and 
approximate  STF's  for  a  static  focus  error  of  0.2  radians  RMS.  For  this  small  focus  error,  the 
exact  STF  is  very  well  approximated  by  the  product  formula  given  by  Eq.  (2.13),  with  the 
approximations  being  numerically  indistinguishable  from  the  exact  STF’s  for  d/ro  =  1  and 
2.  In  Fig.  5b,  the  same  plots  are  shown  for  a  static  focus  error  of  0.7  radians  RMS.  Here,  the 
exact  STF's  are  again  well  approximated  for  d/ro  =  1  and  2  with  a  maximum  relative  error 
of  about  15  percent  at  a  frequency  of  0.6.  For  d/ro  =  4,  however,  the  maximum  relative 
error  is  about  80  per  cent.  In  Fig.  5c,  finally,  the  static  focus  error  has  been  increased  to 
1  radian  RMS.  and  the  approximation  is  no  longer  very  useful.  For  this  case,  the  relative 
error  is  about  40  percent  at  midband  for  d/ro  =  1  and  over  75  percent  for  d/ro  =  2.  For 
d/ro  =  4.  the  approximation  fails  rather  spectacularly  beyond  a  spatial  frequency  of  about 
0.3.  Unfortunately,  accurate  estimates  of  the  STF  are  likely  to  be  most  useful  exactly  when 
d/ro  is  large. 

The  corresponding  STF  results  for  residual  phase  errors  due  to  anisoplanatism  are  il¬ 
lustrated  in  Fig.’s  6  and  7.  In  Fig.  6  the  static  focus  error  is  0.7  radians  RMS,  and  the 
atmospheric  phase  residuals  are  provided  by  angular  offsets  between  the  science  object  and 
guidestar  of  0.5.  1,  and  2  times  the  isplanatic  angle  Bq.  For  B/Oq  =  0.5,  Eq.  (2.13)  provides  a 
good  approximation  to  STF(/c)(,^„P).  In  Fig.  7a,  we  show  STF(k)(o,o)  and  STF(k)(0,,x>); 
along  with  the  individual  factors  STF(K)(o,n)  and  STF(#j)(^„o)  and  the  approximation 
Eq.  (2.13).  Although  STF(K)(o,n)/STF(/c)(o.o)  ~  0.5  at  midband,  the  maximum  relative 
STF  error  using  Eq.  (2.13)  is  still  only  19  percent.  For  the  case  of  B/Bo  =  2  illustrated  in 
Fig.  7b.  however,  Eq.  (2.13)  begins  to  be  unreliable.  The  STF  attenuation  due  to  atmo¬ 
spheric  and  static  sources  are  roughly  equal,  with  STF(k)(,^,,o)  =  0.18  x  STF(<c)(o,o)  and 
STF(/c)(o.i))  =  0.2  X  STF(k)(o.o)  at  midband.  These  conditions  lead  to  a  relative  STF  error 
of  60  percent. 


C.  Gemini  Results 

The  results  obtained  with  the  Gemini-North  telescope  and  AO  parameters  listed  in  the 
final  column  of  Table  1  are  summarized  in  Table  2,  Fig.  8,  and  Fig.  9.  The  formulas  from 
section  2  used  for  these  calculations  are  as  listed  in  subsection  3B  above.  Table  2  list  Strehl 
ratios  at  wavelengths  of  1.20,  1.65,  and  2.20 /im  computed  for  residual  phase  errors  due  to 
atmospheric  turbulence  alone,  optical  aberrations  alone,  and  a  combination  of  the  two  error 
sources.  The  relative  error  in  the  multiplicative  approximation  to  the  overall  Strehl  ratio  is 
no  larger  than  about  four  per  cent  over  this  range  of  cases.  One-dimensional  slices  through 
the  corresponding  OTF’s  for  the  1.65  //m  case  are  plotted  in  Fig.  9.  Once  again,  there  is  very 
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good  agreement  between  the  exact  calculation  of  the  end-to-end  OTF  and  a  multiplicative 
approximation  using  two  separate  transfer  functions. 

In  Fig.  9  we  plot  the  various  STF’s  for  the  “typical”  Gemini  aberrations  and  nominal 
atmospheric  phase  distortions  at  a  wavelength  of  1.2  //m.  The  RMS  phase  is  0.62  radians  for 
the  static  aberrations  and  0.47  radians  for  the  atmospheric  residual.  Because  the  Gemini 
errors  are  not  azimuthally  symmetric,  like  focus,  the  curves  shown  are  radial  averages  of 
two  dimensional  STF's.  With  a  total  Strehl  ratio  of  about  0.3,  the  resulting  PSF  is  within 
the  “sharp  core”  regime,  and  in  practice  speckle  imaging  methods  would  probably  not  be 
required.  Still,  we  point  out  that  since  the  STF  attenuation  due  to  static  aberrations  is  fairly 
evenly  distributed  over  a  range  of  frequencies,  a  relative  error  of  about  8  percent  is  achieved 
at  midband  using  Eq.  (2.13).  This  level  of  agreement  is  obtained  for  a  static  aberration  with 
a  RMS  value  approximately  equal  to  the  fixed  focus  error  considered  in  Fig.  6b,  for  which 
the  accuracy  of  the  multiplicative  approximation  was  much  poorer.  In  comparison  with  the 
OTF  results  presented  in  Fig.'s  4  and  8,  it  appears  that  the  accuracy  of  Eq.  (2.13)  for  the 
STF  depends  more  strongly  upon  the  specific  nature  of  the  static  aberrations. 


4.  Sky  Coverage  Calculations 

The  so-called  “sky  coverage”  of  an  AO  system  for  an  astronomical  telescope  can  be  expressed 
in  terms  of  the  function  P{S  >  So)-  the  fraction  of  a  specified  region  of  the  sky  for  which  the 
system  achieves  an  end-to-end  performance  level  S  at  least  equal  to  a  desired  threshold  Sq.'^ 
Sky  coverage  can  be  computed  for  a  variety  of  performance  metrics,  such  as  time  averaged 
Strehl  ratio  or  the  slit  power  coupling  for  a  spectrometer.®’*®  This  function  is  determined 
by  a  combination  of  guide  star,  atmospheric,  science  instrument,  AO  system,  and  telescope 
parameters,  including  the  unsensed  and  uncorrectable  aberrations  in  the  telescope  optics. 
Sky  coverage  remains  an  issue  for  laser  guide  star  (LGS)  AO  systems,  since  a  natural  guide 
star  (NGS)  is  still  required  to  measure  the  full  aperture  tip  and  tilt  of  the  wave-front.  An 
LGS  AO  system  will  consequently  obtain  best  results  for  science  objects  very  near  to  a 
bright  NGS.  while  performance  for  other  objects  will  be  degraded  by  tip/tilt  jitter  arising 
from  a  combination  of  tip/tilt  sensor  noise  and  tilt  anisoplanatism.  The  AO  system  mod¬ 
eling  techniques  used  to  quantify  these  considerations  in  terms  of  the  sky  coverage  function 
P(S  >  So)  are  described  in  previous  papers.®®  together  with  sample  results  for  the  Gemini- 
North  telescope  which  were  computed  assuming  no  uncorrectable  telescope  aberrations.  In 
this  section  we  investigate  how  these  results  are  influenced  by  the  representative  optical 
aberrations  illustrated  in  Section  3  above. 

The  evaluation  of  AO  system  sky  coverage  for  a  particular  observing  scenario  requires 
a  parameterized  family  of  point  spread  functions  PSF(0,i?)  which  define  AO  system  per¬ 
formance  with  a  natural  guide  star  of  magnitude  R  at  an  angular  offset  0  from  the  science 
field.  Together  with  a  guide  star  model  describing  the  density  of  guide  stars  as  a  function 
of  their  magnitude,  this  family  of  PSF's  can  be  used  to  compute  the  probability  of  locating 
a  natural  guide  star  both  bright  enough  and  near  enough  to  the  science  field  to  achieve  a 
given  level  of  performance.  For  a  LGS  AO  system,  these  point  spread  functions  may  be 
computed  using  the  equation 

PSF(9.  R)  =  PSFho  *  PSFtt(^,  R).  (4.1) 
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Here  PSFj^q  is  the  short-exposure  (tilt  removed)  PSF  due  to  residual  higher  order  aberrations 
left  uncorrected  by  the  laser  guide  star  AO  system,  the  asterisk  denotes  the  convolution 
operator,  and  PSFtt  is  a  Gaussian  transfer  function  which  represents  the  additional  blurring 
in  the  long-exposure  PSF  due  to  residual  tip  and  tilt  errors.  The  short-exposure  PSF  is 
independent  of  the  parameters  6  and  R  because  the  laser  guide  star  used  for  higher-order 
correction  has  constant  brightness  and  can  always  be  projected  in  the  direction  of  the  science 
object.  In  terms  of  this  model,  unsensed  and  uncorrectable  optical  aberrations  impact  AO 
system  sky  coverage  by  degrading  the  point  spread  functions  PSF(0.  i?)  in  two  ways: 

•  The  short-exposure  PSF  and  OTF  are  degraded  directly  according  to  Eq.  (2.6)  above, 
and 

•  Blurring  of  the  NGS  images  in  the  tracker  or  Shack-Hartmann  sensor  used  for  tip/tilt 
sensing  increases  the  tilt  measurement  error  due  to  noise  for  a  fixed  NGS  magnitude 
R.  consequently  increasing  the  residual  tilt  jitter  and  broadening  the  function  PSF^^^. 

For  Gemini,  each  science  instrument  is  equipped  with  its  own  “on-instrument,’’  or  colo¬ 
cated.  Shack-Hartmann  sensor  for  tip/tilt  sensing.  The  magnitude  of  the  increased  tip/tilt 
jitter  due  to  blurring  can  be  quantified  by  evaluating  the  OTF  for  each  sensor  subaperture 
using  Eq.  (2.6).  and  then  applying  standard  formulas  for  quadrant  detector  tilt  measurement 
accuracy  as  a  function  of  signal  level  and  the  shape  of  the  PSF.^^’^®  The  same  aberration 
profile  can  be  used  for  the  tip/tilt  sensor  and  the  science  instrument  since  the  two  devices 
are  colocated.  In  the  calculations  to  date  the  increased  tilt  jitter  due  to  blurring  of  the  NGS 
image  has  been  fairly  negligible,  and  nearly  all  of  the  reduction  in  sky  coverage  associated 
with  static  telescope  aberrations  has  been  due  to  their  direct  effect  upon  the  higher  order 
point  spread  function  PSEj^q. 

Fig.’s  10  and  11  illustrate  sample  results  for  the  effect  of  uncorrectable  aberrations  on 
Gemini-North  sky  coverage.  These  results  are  for  an  imaging  wavelength  of  1.65  pm,  the 
LGS  .AO  and  observing  parameters  listed  in  column  4  of  Table  2,  a  guide  star  density  profile, 
corresponding  to  the  galactic  pole,  and  tip/tilt  sensor  radiometry  parameters  as  described 
previously.*^  Fig.  10  quantifies  sky  coverage  in  terms  of  Strehl  ratio,  while  Fig.  11  considers  the 
case  of  power  coupling  through  a  0.1  arcsec  spectrometer  slit.  In  both  cases  the  effect  of  the 
uncorrectable  aberrations  is  nonnegligible,  but  the  relative  reduction  in  system  performance 
is  smaller  for  the  second  case.  For  example,  the  typical  aberration  profile  reduces  the  Strehl 
ratio  for  50  per  cent  sky  coverage  from  about  0.45  to  0.34,  a  relative  reduction  of  about  24 
per  cent.  The  corresponding  values  for  slit  power  coupling  are  about  0.68,  0.59,  and  13  per 
cent.  Evidently,  the  sample  noncommon  path  aberrations  used  in  these  calculations  contain 
a  range  of  spatial  frequencies  which  scatter  some,  but  not  all,  of  the  central  lobe  of  the  PSF 
outside  of  the  aperture  for  a  0.1  arcsec  slit. 

Many  different  sets  of  observing  conditions  and  noncommon  path  aberrations  may  need 
to  be  considered  for  the  design  and  optimization  of  an  AO  system,  and  it  may  become 
computionally  cumbersome  to  repetitively  evaluate  Eq.  (2.6)  for  each  pair.  Computation 
requirement  can  be  reduced  by  applying  Eq.  (2.8)  or  (2.10)  to  approximate  the  effect  of  non¬ 
common  path  aberrations  on  the  end-to-end  Strehl  ratio  or  OTF.  Fig.’s  12  and  13  illustrate 
the  effect  of  these  approximations  on  the  final  sky  coverage  predictions  for  the  AO  system. 
Using  the  multiplicative  OTF  approximation  leads  to  very  modest  errors  in  the  sky  cov¬ 
erage  estimates  which  are  probably  small  relative  to  other  uncertainties  in  the  calculation. 
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The  errors  in  sky  coverage  associated  with  the  multiplicative  Strehl  ratio  approximation  are 
significantly  larger,  particularly  in  Fig.  13  where  the  approximation  is  applied  to  slit  power 
coupling  instead  of  the  Strehl  ratio. 


5.  Summary 

In  this  paper  we  have  investigated  the  interaction  of  static  telescope  aberrations  and  residual 
turbulence-induced  phase  distortions  upon  the  performance  of  adaptive  optical  systems. 
Multiplicative  approximations  for  the  combined  effect  of  these  two  wave- front  error  sources 
upon  the  time  averaged  Strehl  ratio  and  OTF  of  the  AO  system  are  remarkably  accurate 
as  long  as  the  overall  time  average  Strehl  ratio  is  greater  than  about  0.1.  This  result  allows 
the  effects  of  static  telescope  aberrations  to  be  easily  included  in  AO  system  sky  coverage 
calculations.  A  similar  approximation  for  the  speckle  transfer  function  can  be  seriously 
in  error  unless  the  Strehl  ratio  is  greater  than  about  0.3,  at  which  point  speckle  imaging 
methods  are  generally  no  longer  required. 
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TABLES 


Scenario' 

1 

2 

3 

4 

Aperture  diameter,  m 

1.0 

1.0 

1.0 

7.9 

Linear  obscuration  ratio 

0.0 

0.0 

0.0 

0.152 

Imaging  wavelength,  fira 

1.06 

1.06 

1.06 

1.65 

AO  guide  star 

NGS 

NGS 

NGS 

LGS  at  ^  =  92.5  km 

WFS  Subapertures 

'8x8 

8x8 

8x8 

12  X  12 

DM  Actuators 

9x9 

9x9 

9x9 

13  X  13 

AO  loop  bandwidth,  Hz 

oc 

oc 

oc 

100 

Turbulence  Profile 

HV5/7 

HV5/7 

HV5/7 

Typical  Mauna  Kea,  ^  =  0° 

Static  Aberrations 

Focus 

Focus 

Focus 

Simulated  Gemini- 
North  aberrations 

Fitting  error,  d/ro 

1-4 

1 

1 

0.67 

WFS  noise.  ere/{X/d) 

0.0 

0.0-0.4 

0.0 

0.0 

.4nisoplanatism,  0/6o 

0 

0 

0-2 

0 

Table  1.  Parameter  values  for  OTF  and  Strehl  ratio  calculations.  Here  NGS  denotes  a  natural 
guide  star.  LGS  an  artificaUy  generated  laser  guide  star  at  range  z,  ^  is  the  zenith  angle,  d  is 
the  width  of  a  WFS  subaperture,  tq  is  the  turbulence-induced  atmospheric  coherence  diameter, 
ere  is  the  RMS  WFS  subaperture  tilt  measurement  accuracy,  A  is  the  imaging  wavelength,  9  is  the 
angular  offset  between  the  guide  star  and  the  science  object,  and  6o  is  the  isoplanatic  angle.  The 
notation  “HVs/7’'  denotes  the  Hufnagel- Valley  turbulence  profile  scaled  to  yield  tq  =  5cm  and 
^0  =  7  /irad  at  A  =  0.5 //m.  Futher  description  of  these  parameters  is  contained  in  the  text. 
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Phase  Errors 

Computation 

1.25  //m 

Strehl  ratio 
1.65  //m 

2.20  ;/m 

Turbulence 

Exact 

0.621 

0.856 

Typical  aberrations 

Exact 

0.568 

0.818 

Best  case  aberrations 

Exact 

0.762 

0.913 

Turbulence  and  typical  aberrations 

Exact 

Multiplicative 

0.365 

0.353 

0.545 

0.538 

Turbulence  and  best  case  aberrations 

Exact 

Multiplicative 

0.481 

0.473 

0.651 

0.648 

0,783 

0.782 

Table  2.  Strehl  ratio  results  for  Gemini-North  telescope  and  AO  system  parameters.  This  table  Usts 
Strehl  ratios  computed  in  three  astronomical  imaging  bands  based  on  the  Gemini  system  parame¬ 
ters  listed  in  the  last  column  of  Table  1.  The  first  section  of  the  table  lists  Strehl  ratios  computed 
separately  for  the  effects  of  atmospheric  turbulence  and  optical  system  aberrations.  The  remaining 
two  sections  of  the  table  lists  Strehl  ratios  for  the  combined  effect  of  these  two  error  sources,  where 
the  “exact*"  values  have  been  computed  using  Eq.  (2.7),  and  the  “multiplicative”  values  are  derived 
from  the  first  section  of  the  table  using  Eq.  (2.9).  The  multiplicative  approximation  is  accurate  to 
within  a  few  per  cent  over  the  range  of  cases  considered. 


FIGURES 


Fig.  1.  Representative  unsensed  and  uncorrectable  mirror  figure  errors  for  the  Gemini-North 
Telescope.  This  phase  error  profile  yields  a  Strehl  ratios  of  about  0.66  at  a  wavelength  of  1.65  fxm  , 
which  corresponds  roughly  to  the  Gemini  specification  for  mirror  figure  errors  and  misalignments. 
The  profile  is  a  sum  of  contributions  due  to  primary  and  secondary  mirror  polishing  and  mounting 
errors,  together  with  non  common  path  aberrations  between  the  WFS  and  the  science  instrument. 
The  mirror  figure  errors  have  been  studied  in  considerable  detail  by  Gemini,  but  the  noncommon 
path  errors  used  here  are  merely  “typical”  values  selected  to  match  the  overall  Strehl  ratio  given 
above.  A  more  benign  profile  corresponding  to  an  overall  Strehl  ratio  of  about  0.82  has  also  been 
considered. 

Fig.  2.  Strehl  ratio  calculations  for  representative  AO  system  parameters  and  wave-front  errors. 
These  curves  plot  the  exact  Strehl  ratio  and  its  multiplicative  approximation  computed  for  the  AO 
system  parameters  and  static  and  varying  wave-front  errors  summarized  in  columns  1  through  3 
of  Table  1.  Subfigure  (a)  plots  results  for  wave-front  errors  dominated  by  DM  fitting  error  with 
d/ro  =  1  (top).  2,  and  4  (bottom),  subfigure  (b)  is  for  the  case  of  WFS  noise  with  (Tgf{Xfd)  =  0 
(top).  0.1.  0.2,  and  0.4  (bottom),  and  subfigure  (c)  is  for  anisoplanatic  wave-front  errors  with 
O/Sq  =  0  (top).  0.5.  1,  and  2  (bottom).  In  subfigure  (a),  the  results  for  the  multiplicative  approxi¬ 
mation  overlay  the  exact  Strehl  ratios  to  within  the  resolution  of  the  plot.  See  the  caption  to  Table 
1  for  the  definitions  of  these  variables. 


Fig.  3.  OTF  reduction  due  to  a  static  focus  aberration  of  0.7  radians  RMS. 

Fig.  4.  OTF  calculations  for  representative  AO  system  parameters  and  wave-front  errors.  These 
OTF's  correspond  to  the  AO  system  parameters  and  turbulence-induced  wave-front  errors  de¬ 
scribed  in  the  caption  to  Fig.  2  above,  combined  with  a  static  focus  error  with  a  magnitude  of  0.7 
radians  RMS. 

Fig.  5.  STF  calculations  for  wave- front  fitting  error  with  d/ro  =  U  2,  and  4,  and  the  AO  system 
parameters  listed  in  column  1  of  Table  1.  Exact  STF’s  and  multiplicative  approximations  are 
shown  for  a  static  focus  error  of  0.2  radians  RMS  (Fig.  5a),  0.7  radians  RMS  (Fig.  5b),  and  1 
radian  RMS  (Fig.  5c). 


Fig.  6.  STF  calculations  for  anisoplanatic  wave-front  errors  with  9/6o  =  0.5,  1,  and  2,  and  the 
AO  system  parameters  listed  in  column  2  of  Table  1. 

Fig.  7.  STF's  for  diffraction-limited  imaging,  residual  atmospheric  wave-front  error  only,  and 
static  focus  wave-front  error  only  cases.  The  atmospheric  wave-front  error  is  due  to  anisoplanatism 
with  0/00  =  1  (Fig.  7a)  and  0/0o  =  2  (Fig  7b).  AO  system  parameters  are  as  listed  in  column  3 
of  Table  1.  The  static  error  is  due  to  a  0.7-radian  RMS  focus  error.  Exact  and  approximate  STF’s 
are  also  shown  for  the  combination  of  the  two  effects. 
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Fig.  8.  OTF  calculations  for  Gemini-North  AO  system  parameters  and  wave-front  errors.  These 
OTF's  correspond  to  the  Gemini  AO  system  parameters  and  turbulence-induced  phase  errors 
summarized  in  the  final  column  of  Table  1,  combined  with  typical  and  best-case  unsensed  and 
uncorrectable  optical  system  aberrations.  The  wavelength  is  1.65 /im.  The  OTF’s  are  very  well 
approximated  by  the  product  of  terms  computed  separately  for  the  two  error  sources. 

Fig.  9.  STF’s  for  diffraction-limited  imaging,  residual  atmospheric  wave-front  error  only,  and 
static  telescope  wave-front  error  only  cases.  The  atmospheric  wave-front  error  is  due  to  nominal 
residual  phase  at  1.2fxxa  for  the  Gemini  AO  system  parameters.  The  static  error  is  due  to  the 
“typical  case”  Gemini  optical  aberrations.  Exact  and  approximate  STF’s  are  also  shown  for  the 
combination  of  the  two  effects.  Because  the  Gemini  aberrations  are  not  symmetric,  all  plots  are 
radial  averages  of  two-dimensional  distributions. 


Fig.  10.  Impact  of  uncorrectable  telescope  aberrations  on  Strehl  ratio  sky  coverage  for  Gem¬ 
ini-North  at  an  observing  wavelength  of  1.65 /zm.  Here  “sky  coverage”  is  defined  as  the  fraction  of 
the  sky  close  enough  to  a  sufficiently  bright  guide  star  for  the  AO  system  to  deliver  the  indicated 
Strehl  ratio.  The  uncorrectable  telescope  aberrations  used  for  these  calculations  are  as  illustrated 
in  Fig.  1,  and  the  AO  system  and  observing  parameters  are  listed  in  column  4  of  Table  2.  Infor¬ 
mation  on  the  guide  star  density  model  and  tip/tilt  sensor  radiometry  parameters  can  be  found  in 
a  previous  paper.® 

Fig.  11.  Impact  of  uncorrectable  telescope  aberrations  on  slit  power  coupling  sky  coverage  for 
Gemini-North.  This  figure  is  similar  to  Fig.  10,  except  that  AO  system  performance  is  described 
in  terms  of  the  fraction  of  energy  from  a  point  source  which  is  transmitted  through  a  0.1  arcsec 
spectrometer  slit. 


Fig.  12.  Effect  of  approximations  on  estimated  Strehl  ratio  sky  coverage  for  Gemini-North. 
This  figure  illustrates  how  sky  coverage  predictions  for  Gemini-North  are  effected  if  the  end-to-end 
performance  of  the  AO  system  is  estimated  using  one  of  the  two  approximations  given  at  the  end  of 
Section  2  above.  The  observing,  AO  system,  and  guide  star  parameters  used  for  these  calculations 
are  the  same  as  for  Fig.  9  above,  along  with  the  larger  of  the  two  values  for  uncorrectable  telescope 
aberrations. 

Fig.  13.  Effect  of  approximations  on  estimated  slit  power  coupling  sky  coverage  for  Gemini-North. 
This  figure  is  similar  to  Fig.  12,  except  that  AO  system  performance  is  described  in  terms  of  the 
fraction  of  energy  from  a  point  source  which  is  transmitted  through  a  0.1  arcsec  spectrometer  slit. 
The  results  labelled  “multiplicative  Strehl  ratio  approximation”  were  obtained  by  multiplying  the 
slit  power  coupling  without  telescope  aberrations  by  the  Strehl  ratio  due  to  these  aberrations  alone. 
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Radians  of  phase  at  0.5  p  m 


Figure  1 
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Strehl  ratio 
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Exact  and  approximate  Strehl  ratio  with  fitting  error 


Figure  2a 


Strehl  ratio 


Exact  and  approximate  Strehl  ratio  with  WFS  noise 


Figure  2b 
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Strehl  ratio 


Exact  and  approximate  Strehl  ratio  with  anisoplanatism 


Figure  2c 
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Modulus  of  mean  OTF 


Figure  4a 


23 


Modulus  of  mean  OTF 


Figure  4b 
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Modulus  of  mean  OTF 


Figure  4c 
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speckle  transfer  function 


speckle  transfer  function 


speckle  transfer  function 
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Exact  and  approximate  OTF’s  for  Gemini  scenarios 


speckle  transfer  function 


Fractional  sky  coverage 


Figure  10 


Fractional  sky  coverage 


Figure  12 
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Fractional  sky  coverage 


Figure  13 
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The  work  described  in  this  document  was  accomplished  by  Mr  David  Tyler  in  collaboration 
with  Dr  Brent  Ellerbroek  of  the  Air  Force  Research  Laboratory  (AFRL)  Starfire  Opti¬ 
cal  Range.  Dr  Ellerbroek  funded  Schafer  using  a  continuing  Air  Force  Office  of  Scientific 
Research  (AFOSR)  grant.  The  grant  supports  an  AFRL  collaboration  with  the  Gemini 
Telescopes  Project  {http  :  f lwww.gemini.edu)  to  develop  and  analyze  adaptive  optics  (AO) 
methods,  instrumentation,  and  performance  for  the  Gemini-North  telescope  on  Mauna  Kea 
on  the  island  of  Hawaii.  The  scope  of  Schafer's  work  includes  theory  and  analysis,  numer¬ 
ical  analysis  and  code  development,  and  documentation  of  results  for  journal  publication. 
This  report  summarizes  work  accomplished  between  March  1998  and  November  1998,  when 
.Mr  Tyler  left  Schafer  Corporation.  Attached  is  a  draft  journal  article  detailing  the  work 
performed  in  this  period;  the  acknowledgments  contain  a  statement  that  Mr  Tyler  was  em¬ 
ployed  by  Schafer  when  the  work  was  performed.  The  draft  article  has  been  anonymously 
reviewed  by  two  referees,  and  all  issues  raised  by  the  referees  have  been  addressed  in  the 
attached  version.  Following  is  a  brief  summary  of  the  work  performed. 
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Summary:  Modelling  the  combined  effect  of  static  and  varying  phase  distortion 
on  the  performance  of  adaptive  optical  systems.  Adaptive  optics  (AO)  systems  are 
now  in  use  at  several  astronomical  observatories  throughout  the  world  to  compensate  the 
distorting  effect  of  the  Earth's  atmosphere.  Science  results  from  working  AO  systems  at 
the  William  Herschel,  ESO  3.6m,  and  Canada-France-Hawaii  telescopes  are  encouraging 
enough  that  AO  retrofits  to  such  venerable  instruments  as  the  Hale  200  in.  and  Lick  120  in. 
telescopes  are  well  under  way.  However,  the  primary  mirror  diameters  of  these  telescopes 
range  from  3  to  5  meters,  in  what  is  now  called  the  “intermediate”  range  of  telescope  size. 
The  newest  telescopes  such  as  the  Kecks,  Suburu,  and  the  Gemini  North  and  South  are 
all  substantially  larger,  and  as  a  consequence,  substantially  more  expensive.  Further,  AO 
compensation  of  turbulence  effects  at  optical  and  near-infrared  wavelengths  becomes  more 
difficult  with  increasing  aperture  size.  To  fully  exploit  these  large  telescopes,  then,  AO 
modelling  and  design  must  be  done  at  a  level  of  detail  previously  unnecessary. 

The  existing  literature  suggests  that  most  adaptive  performance  analysis  is  partitioned 
into  separate  treatments  of  time-varying  atmospheric  aberrations  and  static  (or  quasi-static) 
aberrations  in  the  telescope  optics  themselves.  Overall  system  performance  is  typically  esti¬ 
mated  by  multiplying  the  associated  Strehl  ratios  or  optical  transfer  functions  (OTFs).  For 
.AO-compensated  residual  aberrations,  this  multiplicative  approximation  is  not  formally  ex¬ 
act.  Since  both  the  costs  and  expectations  for  large,  ground-based  telescopes  are  so  challeng¬ 
ing.  it  is  prudent  to  examine  in  detail  the  domain  in  which  the  multiplicative  approximation 
is  actually  valid. 

This  summary,  and  the  attached  paper,  describes  such  a  study.  In  particular,  the  ap¬ 
proximations 
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and 


T’{k)(m)«  T{K)tT{K)s/T{KU  (3) 

were  examined,  where  S  is  the  Strehl  ratio,  H{k)  is  the  OTF,  and  T{k)  is  the  speckle 
transfer  function  (STF).  The  subscipts  t  and  s  indicate  time- varying  and  static  quantities, 
respectively.  The  above  quantities  were  calculated  using  numerical  methods  and  compared 
over  a  range  of  static  aberrations  and  a  variety  of  atmospheric  phase  aberrations. 

For  a  reasonably  robust  sampling  of  residual  atmospheric  and  static  aberrations,  the 
multiplicative  approximation  for  the  total  Strehl  ratio  S(t,s)  is  remarkably  good  as  long  as 
that  quantity  is  larger  than  about  0.1.  At  lower  performance  levels,  the  approximation  can  be 
used  as  a  moderately  conservative  lower  bound.  Similar  conclusions  apply  to  approximating 
a  total  OTF  as  the  product  of  two  factors.  For  higher-order  moments  of  the  OTF  such 
as  the  STF.  the  approximation  begins  to  fail  at  somewhat  higher  Strehl  ratios:  however, 
reasonable  accuracy  can  still  be  achieved  if  the  static  aberration  power  is  distributed  fairly 
evenly  over  a  range  of  spatial  frequencies. 

The  validity  of  these  approximations  simplifies  the  design  and  use  of  AO  systems  for 
large  telescopes  as  follows:  During  system  design,  an  overall  system  OTF  may  be  specified 
by  separate  characterizations  of  the  two  error  sources  (static  and  dynamic).  Also,  when  AO- 
compensated  imagery  is  processed  after  detection  (or  “deconvolved”)  to  improve  resolution 
further,  the  total  OTF  can  be  approximated  by  data  from  the  AO  wavefront  sensor  (see 
Ref.  3  in  the  attached  paper)  and  an  optical  system  OTF  as  calibrated  by  a  local  reference 
source. 
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